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Abstract. The RS CVn binary system HR 1099 was extensively observed by the XMM-Newton observatory in February 
2000 as its first-light target. A total of 570 ks of exposure time was accumulated with the Reflection Grating Spectrometers 
(RGS). The integrated X-ray spectrum between 5-35A is of unprecedented quality and shows numerous features attributed to 
transitions of the elements C, N, O, Ne, Mg, Si, S, Fe, and Ni. We perform an in-depth study of the elemental composition of 



Of 



the average corona of this system, and find that the elemental abundances strongly depend on the first ionisation potential (FIP) 
of the elements. But different from the solar coronal case, we find an inverse FIP effect, i.e., the abundances (relative to oxygen) 
increase with increasing FIP. Possible scenarios, e.g., selective enrichment due to Ne-rich flare-like events, are discussed. 
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1 . Introduction The X-ray regime contains a rich forest of atomic spectral 

lines emitted by many chemical species at ionisation stages cor- 

The outer atmospheres of stars are of fundamental importance reS p 0nd i ng t0 p i asma temperatures between ~ 1 and 30 Million 

to the chemical enrichment of interstellar space through stel- K So far> hoW ever, atte mpts to determine the composition of 

lar winds and the ejection of energetic particles (Meyer 1985). stellar coronae other than the Sun > s from x . ray obserV ations 

Interestingly, the elemental composition of the hot magnetic so- (Antunes et al 1994> white 19% Drake 1996) have produced 

lar corona and the solar wind are markedly different from the results that m considere d to be ambiguous, because the low 

Sun's surface composition, which suggests that some fraction- spectral resolution avai i ab l e from even the best nondispersive 

ation process in the denser layers of the solar atmosphere must detectors has not permitted the separation of atomic emission 

selectively enrich the corona with particular elements. Various lines ar ising from different elements. Often, subphotospheric 

observations of the solar corona, the solar wind, solar energetic abundances were claimed . Studies of stellar coronal abun . 

particles, and cosmic rays have established that elements with a dances whh the EIJVE satelHte haye revealed solar . like FIP ef _ 

first ionisation potential (FIP) below 10 eV are overabundant by fects in SQm ^ but nQt aU of me stars studied (Drake et al 1995; 

factors of up to 4 with respect to the solar photosphere, unlike 199?) The high-resolution grating spectrometers on board the 

elements with a high FIP that maintain their photospheric abun- recently launche d American Chandra and European XMM- 

dance (Feldman 1992, see also Meyer 1985). This is usually re- Newton satellites now provide a new opportunity to unambigu- 

ferred to as the "FIP effect". However, the physics of this frac- ously determine both the elemental abundances and the ther- 

tionation mechanism is poorly understood. Comparative abun- mal structure of ste llar coronae through explicit measurement 

dance studies of stars other than the Sun may shed light on the of emission i ine an d continuum intensities in the X-ray band, 

underlying physics. are usm g g ra ting spectroscopy on XMM-Newton to per- 
form such analyses for the first time. The target chosen for the 



Send offprint requests to: Bert Brinkman 



2 



A.C. Brinkman et al.: XMM-Newton RGS observations of HR 1099 



1 — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — T 



0.04 



0.03 - 



0.02 - 



FT 0.01 - + 



0.00 



RGS 1st order spectrum of 
HR1099 




Fig. 1. First order spectrum of HR 1099 
as measured by the RGS instruments on 
board XMM-Newton. The most promi- 
nent spectral features are indicated just 
above the spectrum by the correspond- 
ing emitting ion. 
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present observation is HR 1099, one of the X-ray brightest and 
most magnetically active stellar systems, which was observed 
early in the programme. HR 1099 has been the target of previ- 
ous EUVE and low-resolution X-ray spectroscopic investiga- 
tions (e.g. Pasquini et al. 1989, Griffith & Jordan 1998). Ayres 
et al. (2001) qualitatively compare the Chandra HETG spec- 
tra of HR 1099 and Capella, the latter being a somewhat cooler 
and less active coronal source. A companion paper (Audard et 
al. 2001) to the present work discusses spectral variability in 
HR 1099, based on selected XMM-Newton RGS data sets. 

2. Observations and Analysis 

The XMM-Newton X-ray observatory incorporates a pay- 
load with two identical high-resolution Reflection Grating 
Spectrometers (RGS; den Herder et al. 2001) with a spectral 
resolution of about 0.06 A FWHM. The spectral band cov- 
ered by the RGS ranges from 5 to 35 A. This range contains 
strong lines of C, N, O, Ne, Mg, Si, S, Ca, Fe, and Ni. The 
clearly detected eight consecutive charge states of Fe are of 
crucial importance for the present analysis. HR 1099 was ob- 
served in January and February 2000 for a total of 570ksec 
yielding approximately 1.3 million source counts. The spec- 
trum was extracted by counting events inside a narrow spatial 
window along the dispersed CCD image, and then by using the 
CCD energy resolution to separate the spectral orders. The ex- 
tracted total first-order spectrum between 7 and 35 A is shown 
in Fig. 1. 

3. Results 

In optically thin coronal plasmas, the power of a particular 
spectral line depends on the collisional and radiative transition 
rates. In this work, these quantities are calculated by means of 
the relativistic Hebrew University Lawrence Livermore Atomic 
Code (Bar-Shalom et al. 1998). Owing to the high spectral 



resolution, we can directly scale the contribution of each ion 
species to match the observation. The atomic-level models 
are each calculated for a single electron temperature, T max , 
at which the relevant ion has its maximum abundance. Fig. 2 
shows the theoretical spectrum from 7-18 A (red curve), which 
also includes a phenomenological continuum component (blue) 
corresponding to a bremsstrahlung spectrum (which dominates 
over any other continuum components at the temperatures in 
question), compared with the HR 1099 data points (black). As 
can be seen, the agreement is very good. The majority of lines 
in this region are due to highly ionised Fe as demonstrated by 
the individual ion curves at the bottom of the figure. The con- 
tribution of each ion species yields a direct measurement of the 
amount of coronal gas in the conditions appropriate for that ion. 

With the distance d (= 29 pc) to HR 1099, the calculated 
line power and the measured line flux can be used to ob- 
tain the emission measure EM (in cm -3 ), which is a product 
of the electron density and the H-ion density, integrated over 
the emitting volume. We assume that each ion emits mostly 
around the temperature T max . This leads to a first approxi- 
mation of the emission measure distribution. This distribution 
is sufficiently accurate for our purposes if no steep gradients 
are present. Finally, in order to plot the EM, we use the ionic 
abundances from ionisation balance calculations (Arnaud & 
Raymond 1992, Mazotta 1998) and assume, as a first step, ele- 
mental abundances of the solar photosphere. The resulting EM 
as a function of the electron temperature T c is shown in Fig. 3 
for ions of C, N, O, Ne, Mg, S, Fe, and Ni. These ions sample 
the HR 1099 hot coronae in the T range from kT = 100 eV 
to 1500 eV and illustrate the monotonic increase of the EM in 
this range. The errors on the temperature for each ion represent 
the range of T c for which the ion fraction exceeds 50% of its 
maximum value, except for the He-like ions 6+ , Ne 8+ , and 
Mg 10+ , for which line ratio methods (Gabriel & Jordan 1969) 
are used to determine T c more precisely and the bars reflect a 
25% uncertainty. 
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Fig. 2. Observed HR 1099 flux spec- 
trum (black dots) compared with the 
total calculated spectrum (red curve). 
The inferred thermal bremsstrahlung 
continuum is represented by the blue 
curve. The eight separate plots at the 
bottom show the relative individual-ion 
contributions of Fe 16+ to Fe 23+ (not to 
scale). The non-Fe lines (O, Ne, Mg, 
and Ni) are indicated explicitly. 



Table 1. Abundance estimates of the HR 1099 coronae, Solar photosphere 6 and Solar corona c relative to O 



Element 


FIP 

(eV) 


HR 1099 
coronae 
(present) 


Solar 
photo- 
sphere 


Solar 
corona 


Solar 
corona/Solar 
photosphere 


HR 1099 
coronae/Solar 
photosphere 


Ne 


21.56 


0.57 


0.15 


0.15 


1.0 


3.8 


N 


14.53 


0.17 


0.12 


0.13 


1.08 


1.4 


O 


13.61 


1.0 


1.0 


1.0 


1.0 


1.0 


C 


11.26 


0.28 


0.47 


0.50 


1.06 


0.59 


s 


10.36 


0.0085 


0.019 


0.024 


1.26 


0.45 


Fc 


7.87 


0.0094 


0.038 


0.16 


4.21 


0.25 


Mg 


7.64 


0.043 


0.045 


0.18 


4.0 


0.95 


Ni 


7.63 


0.00070 


0.0021 


0.0089 


4.24 


0.33 


a Error estimates are 20% for all numbers. 
b Feldman(1992). 



c Feldmanet al. (1992). 



Looking only at the eight Fe data points in Fig. 3, the EM 
distribution appears to be approximately a smooth power law 
over the relevant temperature range as indicated by the solid 
line. The use of Fe alone provides a representation of the EM 
behaviour independent of elemental abundance. Note that the 
power-law EM distribution obtained for Fe is strikingly paral- 
lel to the EM trends suggested by the data points for O, Ne, 
and Mg. However, the distribution is vertically offset for these 
other elements. The corrections needed in order to settle these 
offsets allow us to estimate the abundances relative to one an- 
other. The high degree of overlap in the T e ranges of formation, 
especially for the Ne, Mg, Fe, and Ni ions, makes the sense of 
the relative abundance determinations unambiguous. Table 1 
presents our derived approximate abundances with the solar 
photospheric and solar coronal abundances. The latter should 
be considered as average abundances in various solar coronal 
features in which the individual abundances may vary some- 
what (Henoux 1995). Similarly, our HR 1099 results are aver- 
age coronal abundances during our observations. We note that 
we give abundances relative to the well-determined, abundant, 



and well-behaved (on the Sun) high-FIP element O. Relative 
abundances are appropriate for the study of possible dependen- 
cies on the FIR They also avoid complications with the de- 
tailed modeling of the continuum. We plot in Fig. 4 the coronal 
enrichment (relative to solar photospheric abundances) for the 
HR 1099 coronae as a function of the FIR 

We note in passing that the precise composition of the 
HR 1099 photosphere (except for a recent measurement of the 
Fe abundance of both HR 1099 components, Randich et al. 
1994) is unknown. It is expected that overall, the photospheres 
of nearby field stars such as HR 1099 have a mixture of ele- 
ments similar to the Sun. 

4. Discussion 

The Ne/O ratio found here for the HR 1099 coronae is unusu- 
ally high compared with the solar photospheric value, while the 
N/O ratio is also enhanced and the Mg/O ratio is close to solar 
photospheric. The abundance ratios for the remaining elements 
are well below their solar photospheric values. Compared with 
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Fig. 3. The empirical emission measure 
for: O e +, 7+ , Ne 8+ , Ne 9+ , Mg 10+ , 
Mg 11+ , Fe 16+ through Fe 23+ , Ni 18 +, 
C 5+ , N 6+ and S 13+ (in order from left 
to right for each element). The straight 
line is the best-fit power-law for the Fe 
ions alone. 
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Fig. 4. Enrichment of elemental abun- 
dances in the HR 1099 coronae rela- 
tive to (solar) photospheric values, plot- 
ted as a function of the first ionisation 
potential (FIP). All values are given 
error bars of 20%. Note the increase 
with FIP, exactly opposite the trend ob- 
served in the average solar corona. 



the respective ratios of the solar corona, the contrast is even 
sharper. In fact, HR 1099 reveals systematically a reversed FIP 
effect, which is illustrated in Fig. 4; the coronal enrichment in- 
creases monotonically with the FIP. The only exception may 
be Mg (low FIP). Ne is strongly enhanced, by a factor of 3.8 
compared to the solar photosphere. This enhancement is rem- 
iniscent of a similar anomaly observed in a subset of solar 
flares (Schmelz 1993, Murphy et al. 1991). Interestingly, flares 
are thought to be of great importance to the heating of ap- 
parently quiet coronae of active stars (Doyle & Butler 1985, 
Giidel 1994, Audard et al. 2000) and have been suggested to be 
pivotal for the elemental fractionation process as well (Wang 



1996, Giidel et al. 1999, Mewe et al. 1997). Continuous mix- 
ing of photospheric and coronal material in low-lying loops 
through frequent flares could suppress the solar-like FIP effect 
(Feldman & Widing 1990) in these stars altogether. 

It is generally believed that the fractionation is the result 
of diffusive processes across magnetic field lines somewhere 
in the chromosphere, where the ionisation fraction differs sub- 
stantially among elements of high or low FIP (von Steiger & 
Geiss 1989). Low FIP elements can then be accelerated into 
the quiescent corona, leaving an excess of high FIP elements in 
the chromosphere (Schmelz 1993). If this chromospheric layer 
can be accessed by flares (Schmelz 1993), then a high-FIP- 
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rich plasma may be heated and lifted into dense coronal X-ray 
loops giving rise to the reversed coronal FIP effect observed in 
HR 1099. In any case, it is clear that the simple mechanism in- 
volving steady-state diffusion, which has been proposed for the 
quiet solar corona, does not operate in more active stars such as 
HR 1099. Since cosmic rays do show a low FIP enhancement, 
it is unlikely that such active stars can be major contributors to 
the cosmic ray flux. The wide survey of stellar coronae planned 
with the RGS on XMM-Newton will provide more general an- 
swers to the questions raised by these observations. 
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